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ABSTRACT 


This  report  addresses  the  problem  of  estimating  the 
spectral  parameters  of  an  observed  doppler  velocity  spectrum 
using  a  burst  radar  waveform  of  arbitrary  lergth.  Maximum- 
likelihood  theory  is  applied  and  the  exact  M.j.E.  algorithm  for 
estimating  the  spectral  mean  is  derived.  This  M.L.E.  algorithm 
is  shown  to  include,  as  a  special  case,  the  spectral  mean  esti¬ 
mator  originally  proposed  by  R.  W.  Miller  [ 1 ] for  processing  of 

( 

burst  waveforms.  Also,  when  the  burst  waveform  is  a  simple 
pulse  pair,  the  M.L.E.  algorithm  reduces  to  the  spectral  mean 
estimator  originally  proposed  by  W.  D.  Rummler  (2] .  The 
Cramer-Rao  bound  for  estimating  the  spectral  mean  using  burst 
waveforms  is  also  derived.  Simplifications  to  the  exact  maxi¬ 
mum-likelihood  algorithm  are  proposed  and  the  performance  of 
various  estimators  is  compared  to  the  Cramer-Rao  lower  bound. 
Some  preliminary  results  of  studies  of  spectral  width  estima- 
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I. 


INTRODUCTION 


Estimation  of  the  spectral  moments  of  observed  signal 
waveforms  is  currently  an  important  topic  in  many  fields  of 
research.  For  example,  the  analysis  of  pulse  doppler  radar  or 
sonar  echoes  from  distributed  scatterer  targets  is  an  important 
area  of  interest.  Early  work  in  this  area  was  done  by  Rummler 
[2]  who  proposed  estimators  of  the  spectral  mean  and  variance 
from  pulse  pair  radar  waveforms.  His  algorithms  were  based  on 
covariance  argument  techniques  and  are  generally  referred  to  as 
time-domain  processing.  Rummler' s  proposed  algorithms  were 
subsequently  shown  to  be  maximum-likelihood  Estimators  of  spec¬ 
tral  mean  and  variance  (see  K.  Miller  and  M.  Rochwarger f 3 ] ) . 
Processing  of  independent  pulse  pairs  requires  a  long  radar 
observation  time,  and  thus,  it  was  suggested  by  Berger  and 
Groginsky  [4]  that  time-domain  estimation  of  spectral  moments 
fro:?  adjacent  pulse  pairs  of  a  radar  burst  waveform  is  the  most 
direct  and  efficient  technique  to  implement  on  conventional 
radars.  Since  the  explicit  solution  to  the  maximum-likelihood 
equation  was  not  available,  processing  of  the  adjacent  pulses 
of  the  contiguous-pulse  burst  waveform  was  done  simply  using 
Rummler' s  proposed  time-domain  algorithms.  Estimates  from  con¬ 
tiguous  pulse  pairs  with  a  common  pulse  between  each  two  pairs 
were  analyzed  in  Benham  [5]  where  the  performance  accuracy 
(i.e.,  standard  deviation  of  the  spectral  mean  and  variance 


estimators)  was  determined  and  compared  with  independent  pulse 
pair  processing.  More  recent  studies  of  spectral  moment  esti- 
mators  using  contiguous  correlated  pulse  pairs  are  presented  in 
Zrnic  (6,7).  These  references  provide  an  excellent  survey  of 
the  topic  of  spectral  moment  estimation  using  time-domain  pro¬ 
cessing  and  present  a  fairly  complete  mathematical  approach  for 
determining  the  theoretical  accuracy  of  covariance  argument 
estimators. 

It  was  recognized  by  R.  Lee  et  al .  [8]  and  by  R. 
Srivastava  et.  al.  [9]  that  useable  spectral  moment  estimates 
could  be  obtained  by  processing  lags  other  than  the  first-order 
lag  utilized  by  Rummler  and  Berger  and  Groginsky.  Also,  im¬ 
proved  estimates  could  be  obtained  by  averaging  the  individual 
moment  estimates  obtained  at  various  lags.  Lee  proposed  a 
poly-pulse  pair  processing  scheme  for  combining  estimates  ob¬ 
tained  at  various  lags  and  demonstrated  that  by  utilizing 
higher  order  lags  to  obtain  spectral  estimates,  performance  of 
the  estimators  could  be  extended  to  lower  signal-to-noise  ratio 
situations. 

Up  to  this  point  we  have  discussed  spectral  moment 
estimation  using  strictly  time-domain  signal  processing  tech¬ 
niques.  An  alternative  method  of  estimating  the  spectral  mean 
and  variance  utilizing  autocorrelation  estimates  at  all  availa¬ 
ble  lags  combined  with  frequency  domain  processing  was  proposed 


by  R.  W.  Miller  [1].  An  efficient  implementation  of  Miller's 
approach  utilizes  an  FFT  algorithm  and  an  appropriately  select¬ 
ed  window  function.  In  this  alternative  approach,  the  power 
spectral  density  of  the  observed  random  process  is  computed 
from  the  sample  autocorrelation  function  of  the  observed  data 
and  the  spectral  mean  is  then  estimated  to  be  the  location  of 
the  maximum  peak  of  the  computed  power  spectrum. 

In  Section  II  of  the  report,  maximum-likelihood  theory 
is  applied  to  the  problem  of  spectral  moment  estimation  using 
burst  waveforms.  The  exact  M.L.E.  algorithm  for  estimating  the 
spectral  mean  is  derived.  This  optimal  estimator  utilizes 
autocorrelation  estimates  at  all  available  lags  and  may  be 
implemented  either  as  a  time-domain  or  frequency-domain  al¬ 
gorithm.  For  contiguous-pulse  burst  waveforms  the  time-domain 
implementation  of  this  optimal  M.L.E.  estimator  is  shown  to  be 
a  generalization  of  Rummler's  algorithm.  When  the  observed 
data  is  comprised  of  independent  pulse  pairs  the  time-domain 
estimator  reduces  to  the  spectral  mean  estimator  originally 
proposed  by  Rummler  [2] .  The  frequency-domain  implementation 
of  the  M.L.E.  estimator  is  also  shown  to  be  a  generalization  of 
Miller's  algorithm  [1],  In  either  case  (time-domain  or 
frequency-domain  processing)  the  M.L.E.  estimator  utilizes  an 
optimally  weighted  combination  of  the  autocorrelation  estimates 
at  all  possible  lags. 


In  Section  lit  of  the  report  we  present  an  interesting 
analytical  solution  for  the  Cramer-Rao  lower  bound  for  estima¬ 
tion  of  the  spectral  mean  from  pulse  burst  data. 

Section  IV  of  the  report  presents  the  results  of  a 
comparison  of  burst  waveform  processing  versus  pulse-pair  pro¬ 
cessing  using  Rummler  time-domain  processing.  Also,  practical 
simplifications  to  the  M.L.E.  estimator  are  presented  and  stud¬ 
ied,  both  analytically  and  through  the  use  of  a  Monte-Carlo 
simulation.  Finally,  preliminary  studies  of  several  algorithms 
for  estimating  the  spectral  spread  are  presented. 

Section  V  provides  a  brief  summary  of  the  report  and 
suggests  the  need  for  development  of  a  more  accurate  algorithm 
for  estimating  the  spectral  spread  parameter. 


II 


MAXIMUM  LIKELIHOOD  ESTIMATION  OP  THE  SPECTRAL  MEAN 


The  return  signal  is  Assumed  to  be  a  discrete-time 
sequence  of  (complex)  radar  video  samples  spaced  at  the  radar 
pulse  interval,  T.  The  return  signal  is  comprised  of  a  time- 
wise  correlated  narrow  band  signal  component  and  a  white  noise 
component,  i.e., 

*k  ■  sJte^w,tT+  k»0,1»...L-1  (1) 

Both  and  n^  are  zero  mean  Gaussian  processes  centered  at 
zero  frequency.  The  true  signal  power  is  denoted  by  S,  while 
the  correlation  coefficient,  P(kT),  depends  upon  the  time 
spacing  between  samples  of  the  sequence  as  well  as  the  spectral 
width  of  the  siqnal  component.  Thus,  the  autocorrelation  func¬ 
tion  of  the  process  z^  is 

R(kT)  a  S3(kT)ejWkT  +  N5k  (2) 


where  R(kT)*E{z. *z, }  and  N  is  defined  as  the  white  noise 

i  x+k 

power  per  sample.  Also,  the  term  6^  is  defined  as  follows 


6 


k 


=  1 


k  =  0 


(3) 


«  0  otherwise 

The  siqnal  power  spectrum  is  assumed  to  be  symmetric 
about  the  spectral  mean  " w"  and  is  Gaussian  in  shape.  Thus, 
the  correlation  coefficient,  P(kT),  is  also  Gaussia:  and  of  the 
form 


where  the  spectral  width  is  denoted  *W*  in  the  above  expres¬ 
sion. 

In  this  section  of  the  report  we  shall  focus  mainly  on 
estimation  of  the  spectral  mean  of  the  signal  component  and 
shall  derive  the  maximum- likelihood  estimate  of  * <•»*  based  upon 
the  sequence  Uo»  *1»  *2» •  •  •  f*L-1  1*  A  brief  summary  of 
the  mathematical  solution  to  the  proposed  estimation  problem  is 
given  in  the  following  paragraphs. 

Given  the  sequence  of  radar  measurements 


zL-i  *  sL-l  e5“(I,'')T  *  nL-i 
we  write  this  data  set  in  matrix  form 


z  »  Ds  +  n 


where  the  diagonal  matrix,  D,  is  given  by  the  following 


w  A  A  ^  V*  1  r’T'S’T’T'.i 


-  jw(0)T 


j»(1)T 


jn(2)T 


L  *ej«(L-l)Tj 

Since  all  the  variables  are  assumed  to  be  Gaussian,  we 
consider  the  covariance  matrix  of  the  measurement  vector,  z_ 
which  may  be  expressed  as  follows: 

E{«  1*’ }  -  D(C  ♦  NI)D*  (8) 

In  the  above  expression,  C  »  fi{s»  s*  *  }  is  the  covariance 
matrix  of  the  signal  vector,  s.  Denoting  the  logarithm  of  the 
(complex)  Gaussian  probability  density  function  as  "In  f(aO", 
we  have 

In  f(x)  -  -«*,D(C+HIj~1D*'  s  -  ln|c+Nl|  -  L  In  2it  (9) 

The  maximum  likelihood  estimate  of  parameter  "«*  is  obtained  by 


taking 


Un  f(*)}  =  0 


(10) 


Since  the  term  Jc+Nl|  is  independent  of  the  parameter 
we  need  only  consider  the  first  term  in  the  above 
log- likelihood  expression.  For  purposes  of  simplicity,  we 
shall  present  the  maximum-likelihood  solution  for  bursts  of 
length  L»2,  3,  and  4.  The  general  solution  for  an  arbitrary 
number  of  pulses  will  then  be  presented. 
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In  the  following  derivation  we  shall  make  use  of  the 
elements  of  the  inverse  convariance  matrix  (C+NI)“1  which  is 
written 


( C+NI ) 


Y1 1  Y1 2  *  *  *  *  Y1 L 

Y12  Y22 . : 

-yil . YLL- 


For  tne  burst  of  length  2,  we  obtain 


(11  ) 


a 

"7uj 


{in  f  (z)  }  =  -  [ ( - j T )  y1; 


e--l  «P+ 


(jT)Y 


12 


Z0Z1* 


(12) 


Setting  the  above  equation  to  zero  yields  the  expression  which 
the  maximum-likelihood  estimate  w" must  satisfy.  When  this  is 
done,  we  obtain 


V 


e"j“T=  z 


V 


3  «T 


(13) 


Clearly  the  above  expression  is  equivalent  to  the  Rummler 
algorithm  for  estimating  the  mean  of  the  power  spectral  density 
function,  i.e.. 


/s  1  i  , 

“  -  t  ARG  iR<T>  1  (14) 

/\ 

where  R(T)  *  zQ*zi  is  the  unbiased  M.L.E.  estimate  of  the 
complex  autocorrelation,  R(T). 

Similarily,  for  a  burst  of  length  3,  given  data  samples 
{zq,  z-| ,  Z2  ] ,  we  obtain 


A'\ 


5a 


$1 


TTo  Un  f  ( z )  }  « 


-  [ (“ jT)  Y^z^z^  +  (-jT)ifnz1*z,e:’lrfr 


12  1  2 


+  {-j2T)Y13z0*z2e"32a)T  +  (jT)Y12z0z1*e:l^ 


(  15) 


+  ( jT) y12z iz2*e 


j“T  + 


(j2T)Yi3z0z2*e 


j  2  wT 


Setting  the  above  equation  to  zero  yields  the  following 
expression: 

2  Y  /\  /\ 

fz  *z  +  z  *z  +  11  z  *z  e“j  le-^  <i>T 

Lzq  z1  +  z1  z2+  —  zQ  z2  e  je 


=  t2021*  +  21  z2  +  TTT  *c^**+:l"r]e:l“T 


which  may  be  written  as  follows 

y  A  A  /\ 

rO/m,  ,  13  -j  utf  1  -j  OIT 

[R(T)  +  ■■  R(2T)e  J  Je  J 

'12 


=  [R*(T)+  — H-  R*(  2T)e^  ^Je^  ^ 

t12 

In  the  above  equation,  we  have  the  definitions 


R(T)  -  2  (zQ*z1+  z^Zj) 


R(2T)  =  ( z  Q*z  2 ) 

Prom  the  above  expressions,  it  is  clear  that  the 


(  16) 


(  17) 


(18) 


(19) 


maximum-likelihood  estimate  of  parameter  "  ui"  is  obtained  from 


(20) 


Ki  *.-  *.~  <.T"i-i  VT. *.1  V.  VIW ».'i VI *.1^"  '-~.,''.,'-’J  '  W,.^ V 


*  V  '  ' 

/\  1  t  T13  ^  o»T  i 

«  *  i  ARGiR(T)  +  — ^  R(2T)e  J  } 

1  t12 

/N 

Since  the  above  solution  for  <*>  is  implicit,  one  iterative 

method  of  solving  the  expression  is  suggested  as  follows. 

Given  the  data  samples  {zq,z-\,Z2  }»  obtain  an  initial 

/s  X 

approximation  of  from  R(T),  i.e.,let 

^  ^  ARG  {1r(  T )  }  (21) 

The  improved  estimate  of  to  based  upon  R(T)  and  R(2T)  is  then 
obtained  iteratively  usinq  the  expression 

/X  ,  Y13  /s  -j^T, 

“i  +  1  T  ARG{R(T)  +  y-2-  R(2T)e  1  }  (22) 

1 2 

For  a  burst  of  length  4,  given  data  samples  {zq ,z ^ ,&2 » *3 1 
we  obtain  the  M.L.E.  estimate 

4  Y  .  ^  v  .  ^ 

/\  1  ’*13  ^  -  -i  jut  1 1 4  ^  _-i  2  (oT  i 

«  =  i  ARG  (r( T)  +T-—  R(2T)e  3  l+  -ii-R(3T)e  }  (23) 

1  JT12  T 1 3 


where  we  have  defined 


R(T)  =  $(*„%,  +  ^  2,*  *2  +  22*  z3) 


(24) 


/N  I  *  * 

R(2T)  =  2  ( zg  z2  +  Z1  z3 ) 


R(3T)  *  (zQ  z 3 ) 


(25) 


(26) 


We  note  from  the  above  equations  that  the  estimate  R(T)  is 


S\1 


ft 


% 

tv 


$ 


k  ' 

k*  . 

SO 

t* 

K* 

a 


biased,  since  in  general,  Y23  *  ^12*  In  a  real-time 
implementation  of  the  M.L.E.  estimator  one  might  use  an 
unbiased  estimate  of  R(T)  in  the  solution. 

Finally  we  present  the  solution  for  the  estimate,  'u't 
for  an  arbitrary  burst  length.  Omitting  the  details,  the  solu¬ 
tion  is  given  by  the  expression 


/V 


L-1  ^  -A. 


«  =  f  ARG  l  I  Ck  R(kT)e~j “(k"1 )T} 
1  k=  1  K 


(27) 


where  the  weighting  coefficients,  C are  defined  in  terms  of 
the  elements,  Y^j,  of  the  inverse  covariance  matrix,  (C+NI)~V 
We  note  that,  in  general,  the  M.L.E.  estimator  utilizes  esti¬ 
mated  autocorrelations 
ble  lags  to  determine  since,  in  general,  the  coefficients 


/\  /N  /\ 

R(T),  R( 2T) , . . . . ,  R(L-1)T  at  all  availa- 


Yij  *  0.  The  estimator  is  implemented  as  a  time-domain  algo¬ 
rithm  and  it  is  seen  (by  comparing  Equations  (14)  and  27))  to 
be  a  generalization  of  the  spectral  mean  estimator  proposed  by 
Rummlerf^j.  This  estimator,  however,  is  implicit  in  struc¬ 
ture  and  an  iterative  scheme  must  be  utilized  to  find  the  'w*' 
which  satisfies  Equation  (27).  Thus,  we  are  interested  in  de¬ 
termining  an  explicit  method  for  obtaining  the  estimate  ^  In 
the  following  paragraphs,  we  develop  a  frequency-domain  imple¬ 
mentation  of  the  estimator  which  is  shown  to  provide  an  explic¬ 
it  solution  for  'u'.  This  alternative  implementation  is  a  gener- 
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-V 

i 
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alization  of  the  spectral  mean  estimator  proposed  by  Miller1  . 

/V 

Prom  the  estimates  of  the  autocorrelation  function  {R(0), 

A  A  , 

R(T) , . . . ,R( L-1 )TJ  we  compute  an  estimated  power  spectral  den¬ 


sity,  S(<*>),  as  follows 


(L-1) 

v  /'N 

S(o)  *  l  A.R(kT)e 

k=-(L-l) 


j  wkT 


(28) 


where  the  coefficients  (Afc)  comprise  an  arbitrary  set  of 

/\  /\ 

weights  or  equivalently  a  window  function.  Since  R(-kT)*R*(kT) 


we  may  write 


<L-1>  /v 

S(«)  *  l  A  tR*(kT)fe 
k*1 


-1  okT 


(29) 


A  (LI)  /\  • 

+  AnR(  0)  +  I  A.RtkTJe3 
°  k*1 

Miller  1 1  ]  took  as  the  estimate  of  spectral  mean,  u>,  that 

✓\  /N  , 

frequency  “>  which  maximized  the  function  S(“>).  That  is,  w  is 
estimated  as  that  frequency  corresponding  to  the  peak  of  the 
power  spectral  density.  Thus,  o'  is  found  by  taking 


s<“>  I  -  0 

/W' 


(30) 


Performing  the  indicated  differentiation  one  may  easily  obtain 
the  result 


■SVV\  V.V.VV.V.  A".*/.--'  .'•>  YXA.  V . •  V\- 


EWswwnpwww^iacKiiJC*  <rv «~.  *V  •'J  «V  ^WWWV.’V  ' * 


L-1 


l  CwR(kT)e‘j  W(k_1  )T  e^"1 
k*  1  * 


(31) 


L-1 


I  C.  Rv(kT)e;ju,(k‘1  )T  e^“T 
k»1  K 


where  we  have  assumed  ■  A_k  and  substituted  C*  »  kA^ 
into  the  above.  But  the  above  expression  may  be  shown  to  be 

i 

( 

equivalent  to  the  M.L.E.  estimate 

A 


/\  i 

U>  =  ^  ARG 


Ck$(kT)e"j“(k_1)TJ 


(32) 


/\ 


Therefore,  instead  of  solving  for  u  from  the  implicit  formula 
of  Equation  (32),  we  may  compute  the  power  spectral  density 


j  function,  S(w)  of  Equation  (28)  and  locate  the  peak  of  the 


function  to  determine  w  explicitly. 
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III.  CRAMER- RAO  BOUND  FOR  SPECTRAL  MEAN  ESTIMATES 

Next,  we  wish  to  determine  the  Cramer-Rao  lower  bound 
on  achievable  estimation  accuracy  by  determining  the  Fisher 
information  matrix  for  the  parameter  "  .  Again  we  shall 

present  the  general  solution  by  first  deriving  the  results  for 
short  bursts  of  length  2  and  3.  The  general  solution  for  a 
burst  of  length  L  will  then  be  presented. 

For  the  case  of  L»2  pulses,  we  obtain 

— 2  {in  f(z)}  «  -[(-jT)2Y1220*z1e“j<*T  +  ( jT) 2  Y^z^  *ej  ]  (33) 


Taking  the  "expectation"  of  the  above  yields  the  Fisher 
information  matrix,  which  is  found  to  be 


{in  f(z)}J»  -2T  Y12S  (5(T) 


(34) 


Substituting  Y12  into  the  above  expression,  one  obtains  the 
Cramer-Rao  lower  bound  for  the  variance  of  the  estimate  ^ 


>  1 


(  1+N/S) 
B2(T) 


(35) 


Incidentally,  the  above  expression  is  identical  to  the 
approximate  predicted  accuracy  of  the  Rummler  pulse  pair  method 
for  estimating  the  spectral  mean,  ^(see  Rummler  (2]).  His 
derivation  involved  approximating  the  nonlinear  estimator 
(Equation  (14))  using  perturbation  theory  and  computing  the 
variance  of  the  linearized  estimate,  w. 
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Similarly#  for  a  three  pulse  burst  we  have 


* 

dw 


{in  f(z)}  -  -[(-jT}2Y12*0%1e“j<‘*r  +  (-jT)2Y12s)#*2e’j“T 


+  (-j2T)2Y13*0**2e_j  (jT)2Y12z0*1*ej 


wT 


(36) 


+  ( jT)2Y12z1z2*  ej<rfr  +  (j2T)2Y13z{)z2*e:i2“T] 

Finally#  for  the  three  pulse  burst#  the  Cramer-Rao  lower  bound 
becomes 


°^2> 


1 


(37) 


4TZ[yi2  SB(T)  +  2Yn  S  0(  2T)  ] 

The  Cramer-Rao  bound  for  a  pulse  burst  of  length  L  is  similarly 
determined.  The  derivation  is  briefly  summarized  as  follows. 
For  the  general  L-pulse  burst  we  write  the  log- likelihood 
function  as  follows: 
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I 


s"' 

4 

*: 


+  CONST  (38) 

where  the  remaining  terms  which  are  independent  of  the 
parameter  "<»>"  are  simply  denoted  as  "CONST".  Before  proceeding 
with  the  derivation  we  make  the  following  remarks.  In  the 
above  expression,  which  is  of  the  form 

In  f(z)  »  -  2*  D(C+NI)-1D  z  +  CONST  (39) 

we  have  assumed  the  matrix  inverse  to  be  more  generally  of  the 
form 


(C+NI)"1  =  (Yij)  (40) 

* 

Also,  we  have  combined  the  matrices  D  and  D  with  the  elements 
of  the  matrix  (Y^j).  Hopefully,  this  will  lead  to  the  desired 
solution  in  a  clear  manner.  Computing  the  required  second 
partial  derivative  with  respect  to  the  parameter  "  u>"  and  then 


«.  *  „  *  ^  » 
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factoring  out  the  complex  exponentials  leads  to  the  result 

1  32  ,  ,,  %  *  *  jwT  *  j2uJT  *  j  (L“l)uiT. 

-a  — j  In  f(*}  -  (*0,l!le  '*2*  '*  L-l  ) 

T  3w 


0  hl2  4  Y13  9  Y14 


Y21  ° 


4  Y31  Y32  ° 


i 9  Y41  4  Y42  Y43 


4  Y24  9  Y25 


-  j  2uT 


YL-1  ,Ij 


9YL,L  -3  4  YL,L-2  j Yl,L-1  I  ° 


-j(L-l)u)T 


L’his  expression  is  then  expanded  as  the  quadratic  form 


7  77 


{in  f< 


2  *  .  (m-n) uT 


m*0  n=»0 


(n-m)  y,  z„  *  e 
'  mn  hi  n 


(42) 


Finally,  taking  the  "Expectation"  we  obtain  the  result 

E  5  {in  f(z)  })  *22  (n-m)2Y  B(  ( n-m)T)  (43) 

T  i  \  i  )  ir«0  n»0 


where  we  have  used  the  fact  that 


e{z*z  }  “  S  8((n-m)T)e^ 
m  n 


.  (n-m) wT 


(44) 


Summarizing  our  derivation,  we  have  found  the  following  general 
solution  for  the  Fisher  Information  matrix 
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(45) 


(-  i-j  tin  f(«)l)  *  -  T2  *£’  (i-j)J'rijAil 

\  /  i«0  j«0  13  13 


where  we  define  the  LxL  matrices 

(A.j)  -  (C+NI) 


(46) 


(  Yi;j )  -  (ONI) 

Then  the  C-R  bound  is  given  by  the  following  expression. 


>  — 
“  T2 


M  1  M  1  «% 

I  /  (i“j)  Y.a  A., 

i-0  j-0  13  13 


(47) 


(48) 
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IV 


ANALYSIS  AND  SIMULATION  RESULTS 


First  the  results  of  a  comparison  of  burst  waveform  processing 
versus  pulse-pair  processing  using  Rummler  techniques  is 
presented.  For  this  comparison,  the  total  number  of  available 
pulses  is  considered  to  be  fixed  (specifically  we  take  L-16 
pulses)  and  compare  estimator  performance  for  a  single  burst 
with  the  estimator  performance  achieved  using  8  independent 
pulse-pairs  (again,  processing  a  total  of  16  radar  pulses). 
Figures  1,  2  and  3  show  the  results  of  this  study.  In  each  of 
these  figures  three  performance  curves  are  shown.  First,  the 
fundamental  accuracy  achievable  with  a  16  pulse  burst  waveform 
is  shown,  based  upon  the  C-R  bound  derived  earlier  in  this 
report  (see  Equation  48).  Also,  the  theoretically  predicted 
accuracy  achievable  by  processing  the  16  contiguous  pulses  of 
the  burst  using  the  Rummler  time-domain  algorithm  is  shown. 
Thus,  we  evaluate  the  accuracy  of  the  estimator  (see  Appendix 
3) 

«  »  ^  ARG  {R(T)  }  (49) 


where 


R(T) 


1 

I7-T 


* 


zizi+1 


(50) 


Finally,  the  accuracy  achieved  oy  Rummler  processing  of  8  in¬ 
dependent  pulse  pairs  (hence,  a  total  of  16  transmitted  radar 
pulses)  is  shown.  This  is  considered  to  be  a  fair  comparison 
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in  that  both  approaches  utilize  the  same  total  number  of  radar 
pulses.  Thus,  both  approaches  utilise  the  same  total  transmit¬ 
ted  radar  energy.  The  major  difference,  of  course,  is  that  the 
8  independent  pulse  pairs  use  a  considerably  longer  radar  ob¬ 
servation  time. 

The  curves  of  Figures  1  through  3  show  normalized  esti¬ 
mation  accuracy  versus  normalized  spacing  between  the  pulses  in 
the  burst  or  pulse  pair.  The  normalization  used  in  these 
curves  is  described  as  follows.  The  assumed  signal  correlation 
function  is  given  by  the  expression 


0(T) 


2  2 
e-2**(WTr 


(51) 


where  W  «  2  0^/i.  we  have  defined  the  signal  process  decorrela 

tion  time  to  be  that  time,  Tq,  required  for  the  signal  process 

-1/2  -1/2 
to  decorrelate  to  e  '  m  0.6  in  value.  Setting  $(TQ)  ■  e  ' 

we  find 


$ 

I 

I 


.  » 


(52) 


Therefore,  in  Figures  1,  2  and  3  we  have  evaluated  the  normal¬ 
ly 

ized  estimation  accuracy  ( ST. DEV {v }/ )  versus  normalized  time 
between  pulses  (T/Tq).  Note  that  the  mean  doppler  velocity 
(m/s)  is  related  to  the  spectral  mean  estimate  through  the  for- 

A  A 

mula  V  *  (A/4*)u>.  The  pulse-pairs  are  assumed  to  be  suffi- 
cently  separated  in  time  so  that  statistically  independent 


\' 
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sampling  of  the  signal  process  from  pair  to  pair  is  achieved. 
The  figures  show  performance  for  three  selected  signal-to-noise 
ratios  (-6  dB,  0  dB,  +10  dB)  and  are  valid  for  any  specific 
spectral  spread.  From  these  curves  one  may  draw  the  following 
conclusions. 

1 .  In  each  case  it  is  clear  that  burst  processing 
utilising  a  total  of  15  "lags  of  1”  available  from 
the  16  contiguous  pulses  is  superior  to  pulse-pair 
processing  using  8  Independent  Rummler  pulse-pairs. 
This  is  especially  true  at  low  signal-to-noise 
ratios.  At  high  signal-to-noise  ratios  both  esti¬ 
mators  provide  comparable  performance#  the  burst 
waveform  having  a  small  advantage. 

2.  Comparing  the  performance  of  the  16  pulse  burst 
waveform  with  its  corresponding  theoretical  C-R 
bound  shows  clearly  that  the  16  pulse  burst  wave¬ 
form  with  Rummler  time-domain  processing  of  the  15 
"lags  of  1"  products  achieves  estimation  accuracy 
comparable  with  the  C-R  bound  except  for  low  sig- 
nal-to-noise  ratios.  Thus#  it  appears  that  some 
improvement  in  estimating  the  mean  velocity  is 
possible  for  only  the  low  (S/N)  ratio  case. 

3.  From  the  curves  shown  in  Figures  1,  2  and  3  it  is 
also  seen  that  the  optimal  pulse  spacing  is  given 
by  T/Tq  *  1#  i.e.#  selecting  the  pulse  spacing 

T  *  Tq  minimizes  estimator  errors  for  both  pulse 
pair  waveforms  and  burst  waveforms.  Thus#  to 
specify  an  optimal  interpulse  spacing  T  ■  Tq  one 
must  have  knowledge  of  the  true  spectral  spread, 
ov#  as  indicated  in  Equation  (52). 


In  the  remaining  paragraphs  we  restrict  our  discussion  to  burst 
waveform  processing  and  consider  the  design  of  optimal  burst- 
forms.  The  basic  objective  will  be  that  of  achieving  good  es¬ 
timator  accuracy  using  a  burst  waveform  having  a  small  inter¬ 
pulse  spacing  (T<Tg).  Note  from  the  curves  of  Figures  1 
through  3  that  excellent  performance  can  be  obtained  with  the 
burst  waveform  provided  the  time  between  pulses,  T,  is  properly 
selected.  At  high  signal-to-noise  ratios,  selection  of  a  good 
pulse  spacing  is  not  difficult  since  the  performance  cuLves  are 
quite  flat,  i.e.,  are  very  insensitive  to  the  parameter  T  over 
a  wide  range  of  values  (see  Figure  3).  Even  with  the  signal- 
to-noise  ratio  as  low  as  a  few  dB  this  insensitivity  to  para¬ 
meter  T  is  true.  Thus,  for  any  reasonable  signal-t(  -noise 
ratio  the  spacing  between  pulses  could  for  example  be  selected 
between  0.3TQ  to  1.5  Tq  and  excellent  performance  will  be 
achieved  over  the  equivalent  range  of  the  spectral  spread  para¬ 
meter.  There  is  no  need  to  improve  upon  the  Rummler  burst  es¬ 
timator  performance  over  this  range  of  parameters.  Therefore, 
we  next  consider  the  problem  of  improving  performance  at  low 
signal-to-noise  ratios. 

For  lower  signal-to-noise  ratios  (say  0  dB)  selection 
of  the  interpulse  spacing  begins  to  become  a  more  important 
consideration  (see  FijUre  2).  For  example,  selecting,  T=0.3TQ 
one  observes  that  burst  waveform  performance  begins  to  degrade 


relative  to  the  estimation  accuracy  predicted  by  the  C-R  bound. 
For  very  poor  signal-to-noise  ratios  (say  -6  dB)  selection  of 
the  interpulse  spacing  begins  to  become  very  critical  (see  Fig¬ 
ure  1)  and  with  T  ■  0.3Tq,  burst  estimator  performance  becomes 
very  degraded  relative  to  the  CR  bound.  For  this  reason  we 
have  considered  the  problem  of  improving  burst  estimator  per¬ 
formance  at  these  very  low  signal-to-noise  ratios. 

We  remark  that  in  deriving  the  MLE  spectral  mean  esti¬ 
mator  (see  Equation  27)  it  was  tacitly  assumed  that  the  covari¬ 
ance  matrix  (C+NI)“1  was  exactly  known  a  priori  or  has  been 
accurately  estimated  from  observed  data.  However,  if  the  coef¬ 
ficients  (Yij)  must  be  estimated  from  the  data  then  applying 
the  maximum-liklihood  approach  to  the  simultaneous  estimation 
of  the  parameter  set  {w,  y  ,  y^2,  . ..,  YLL}  requires  the 
solution  to  the  following  equations 

4m  Un  f<!>  1  =  0 

{in  f (z)  }  =  0  (53) 

{in  f (z)  }  =  0 


3 

7? 
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3 

T? 


LL 


The  solution  to  this  set  of  equations  yields  the  joint  maximum 


liklihood  estimates  {w,  }.  This 

the  estimator  of  " to  be  of  the  form 


of  course  provides 
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1  ARG  £  CkR(kT)e'jW(k~1)T 
A  k=1  * 

where  the  C.  are  defined  in  terms  of  the  M.L.E.  estimates  of 

A  K  ^  « 

the  of  the  inverse  covariance  matrix  (C+NI)~  .  For  an 
L-pulse  burst  this  implies  the  simultaneous  solution  to  1  +  0.5 
L ( L*5- 1 )  equations.  Letting  L*16  this  means  one  must  solve  a 
total  of  137  coupled  nonlinear  equations  -  a  formidable  task! 

The  structure  of  the  M.L.E.  spectral  mean  estimator 
suggests  that  using  the  information  contained  in  the  higher 
order  lags  of  the  burst  waveform  one  may  improve  the  accuracy 
in  estimating  the  spectral  mean  at  low  signal-to-noise  ratios 
for  small  interpulse  spacings.  To  this  end,  we  consider  the 
simpler  approach  first  suggested  by  Lee  (see  Reference  [8]). 

In  this  appproach,  estimates  of  the  spectral  mean  are 

A  A 

obtained  using  each  estimated  autocorrelation  {R(T) ,  R(2T), 

A 

R(3T ),...}  using  Rummler’s  time-domain  method.  In  these 
studies,  we  have  considered  the  following  three  spectral  mean 
estimators 

V,  =  ^  ARG  { R  ( T )  } 

V2  *  THJW)  a  arg{r(*2T) } 

^3  *  TvTSf)  -  arg  {r(  3T )  } 


(55) 


where 


L-2 


A  1  V*  * 

R(T>  '  E=T  S  *1+1 

1=0 


1  L“3  * 

R(2T)  -  .£  *i  *i 

*  ■}  * 
R(3T>  *  L=I  £  Zi  Zi 


+2 


+  3 


(56) 


Then  the  estimates  V1 ,  V2  and  V3  are  combined  in  a  linear 
fashion  to  obtain  improved  performance.  To  this  end,  we  have 
considered  spectral  mean  estimators  of  the  form 

v,2  -  =,v,  ♦  c2v2 


V1 23  *  C1V1  +  C2V2  +  c3v3 


(57) 


where  the  coefficients  c^t  c2,  c3  are  selected  either  to  mini- 

A  A 

mize  the  rms  estimation  error  of  V12,  V123  (optimal  weights)  or 
simply  to  average  the  estimates  (uniform  weights). 

Figure  4  shows  the  theoretical  estimation  accuracy  of 

A  A  A 

each  of  the  spectral  mean  estimators  V1 ,  V2  and  V3,  computed 

s\  /\ 

from  estimated  autocorrelations  R(T),  R(2T)  and  R(3T),  respec¬ 
tively.  For  the  specific  case  considered  (S/N  “  -6  dB, 
ov  =  150  m/s  *=.03m,  L=16  pulses)  the  curves  indicate  theoreti¬ 
cal  accuracy  versus  spacing,  T,  between  adjacent  pulses  of  the 
burst  waveform.  If  it  is  desirable  to  utilize  a  quick  burst  of 

pulses  corresponding  to  a  small  spacing  T,  we  see  for  the  data 

/\ 

shown  m  Figure  4  that  the  spectral  mean  estimator,  V3,  pro- 
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vides  best  performance.  For  pulse  spacings  greater  than  8  u 

/\ 

seconds  the  performance  of  estimator  V-  degrades  rapidly.  How- 

/\  /\ 

;:ver,  in  this  region  of  pulse  spacing  both  estimators  and  V2 
provide  good  performance. 

/s 

Given  the  three  estimators  ,  V2  and  V3 ,  each  of  which 
provide  useable  estimates  of  the  spectral  mean,  it  is  re  tuna¬ 
ble  to  consider  obtaining  improved  estimates  using  a  linear 

combination  of  the  individual  estimates.  Note  however  that  the 
/\  /\  /\ 

estimates  V1f  V2,  and  V3  are  correlated  and  the  amount  of  im¬ 
provement  in  estimation  accuracy  achievable  depends  upon  the 
degree  of  correlation. 

Figure  5  shows  the  results  of  a  theoretical  analysis  of 

/N  /\ 

the  performance  of  the  spectral  mean  estimators  V12  and  V1 23 

using  uniform  weights  as  in  Lee  [8]  (for  comparison  the  figure 

/\ 

also  indicates  the  performance  of  estimator  ) .  The  curves  of 
Figure  6  indicate  the  theoretical  performance  of  these  same 
estimators  using  optimal  minimum  variance  weights.  Comparing 
the  curves  of  Figure  5  and  Figure  6  it  is  seen  that  the  simpler 
method  of  averaging  estimates  using  uniform  weights  provides 
significantly  improved  accuracy  relative  to  the  single-lag  es- 
timator  .  This  is  especially  true  for  small  pulse  spacings 
and  at  large  spacings  for  ^23'  and  at  still  larger  spacings 
for  2 ,  and  V1 .  For  this  case,  the  additional  improvement 
achieved  using  optimal  weights  may  not  be  justifiable.  That 


r  v_  ■.*  V  *." 


is,  since  the  computation  of  optimal  weights  requires  knowledge 
of  both  the  S/N  ratio  and  spectral  width  parameters,  implemen- 
tion  of  optimal  weights  may  not  justify  the  computational 
expense. 

Figure  7  presents  the  results  of  a  Monte-Carlo  simula¬ 
tion  of  spectral  mean  estimators^,  and'v'lj^.  Por  the  case 

considered  ( =  150  m/s,  0.03  m,  L*16  pulses  and  T=5  P  sec) 
we  have  evaluated  estimator  accuracy  as  a  function  of  signal- 
to-noise  ratio.  We  have  also  included  both  the  Cramer-Rao 
lower  bound  and  the  theoretically  determined  accuracy  for 
algorithm  [see  Appendix  B] .  From  the  figure  we  observe  that 
the  simulation  of  estimator  'vl  agrees  quite  well  with  accuracy 
predictions  at  high  signal-to-noise  ratios.  Also,  almost  no 
improvement  in  performance  was  obtained  at  high  signal-to-noise 
ratios  using  estimator  V^2  or  Vj2^.  At  low  signal-to-noise 
ratios  the  performance  of  estimator  departs  from  the  theo¬ 
retical  predictions  due  to  approximations  used  in  the  linear¬ 
ized  error  analysis.  As  expected,  however,  at  low  signal-to- 
noise  ratios  we  observe  a  significant  improvement  in  perform¬ 
ance  achieved  using  estimators  and  V^2j.  Finally  we  remark 

that  uniform  weights  were  used  in  this  simulation. 

Figure  8  shows  the  results  of  a  Monte-Carlo  simulation 
of  several  spectral  width  estimators  which  were  evaluated. 

Among  the  algorithms  simulated  are  Rummler's  original  spectral 
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tions.  Secondly,  the  Cramer-Rao  lower  bound  is  observed  to  be 
very  optimistic  and  does  not  provide  a  tight  bound  on  the  per¬ 
formance  of  these  spectral  width  estimators. 


V 


SUMMARY 


This  report  considered  the  problem  of  estimating  the 
spectral  mean  and  spectral  width  parameters  of  an  observed  dop- 
pier  velocity  spectrum  using  a  burst  waveform  of  arbitrary 
length.  The  observed  waveform  was  modeled  as  a  correlated 
Gaussian  random  process  and  the  maximum  likelihood  algorithm 
for  estimating  the  spectral  mean  was  derived.  Also,  the  cor¬ 
responding  Cramer-Rao  lower  bound  for  estimation  of  the  spec¬ 
tral  mean  parameter  was  derived.  The  M.L.E.  spectral-mean 
estimator  was  shown  to  include  both  the  time-domain  spectral- 
mean  estimator  proposed  by  Rummler  [2]  and  the  frequency  domain 
spectral  mean  estimator  proposed  by  Miller  [1).  Simpl ications 
to  the  M.L.E.  spectral  mean  estimator  were  proposed  and 
studied,  both  theoretically  and  by  using  Monte-Carlo  simula¬ 
tions.  It  was  shown  that  by  using  additional  order  lag  estima- 

/\ 

tes  (e.g.,  R(T) ,  R(2T),  R(3T))  that  improved  spectral  mean 
estimates  could  be  obtained  and  this  improvement  in  performance 
could  be  achieved  at  lower  signal-to-noise  ratio  situations. 

Finally,  preliminary  studies  of  several  spectral  width 
estimators  were  presented  based  upon  Monte-Carlo  simulations. 
From  these  simulation  results  it  was  observed  that  performance 
of  the  spectral  width  estimators  was  far  from  the  Cramer-Rao 
lower  bound,  thus,  there  is  need  for  development  of  a  new  im¬ 
proved  algorithm  for  estimating  the  spectral  width  parameter. 
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APPENDIX  A:  COMPLEX  RADAR  SIGNAL  SIMULATOR 


This  appendix  describes  the  technique  used  to  simulate 
complex  radar  signal  samples  for  burst  waveforms  of  arbitrary 
length.  The  radar  signals  being  simulated  are  described 
earlier  in  the  report  (see  Equations  1  through  6).  The 
technique  has  been  used  to  simulate  burst  waveforms  with  as 
many  as  16  pulses  and  as  few  as  2  pulses.  The  technique  will 
be  described  using,  as  an  example,  a  three  pulse  burst.  The 
extension  to  the  general  L-pulse  case  is  then  trivial. 

Assume  we  have  three  complex  pulse  returns 


S0e 


s.j  e 


s2e 


jo  UT 

+  n 

j  1  iUT 

+  n 

j2uT 

■f  n 

(A-1) 


From  the  assumed  signal  autocorrelation  function,  which  is 

(A-2) 


e{z*  zi+kl  =  S0(kT)ejkU>T  +  N  «k 


we  obtain  the  covariance  matrix  for  the  vector  ss'=(  Zq,z^  ,22) 


z 

+ 

w 

1 _ 

S  6(T)e”^ 

SB(2T)e"j2aiT 

E 1 z  z* '  }  = 

SB(T)e^  “T 

S  +  N 

SPtTle"^ 

S0(2T}e^(1)T 

S  0(T)e^ MT 

S  +  N 

(A-3)  _ 
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The  basic  idea  is  to  recognize  that  this  covariance  matrix  can 
be  factored  into  the  product  of  three  matrices,  i.e. 

Elz  z*'}  «  D(C+NI)D*'  (Pi- A  ) 


1 

0 

0 

Where  D  * 

0 

ej“T 

0 

-  0 

0 

ej2“T  J 

(A-5) 


S  +  N 

S  B(T) 

SB(2T) 

and  ( C4NI  )=. 

SB(T) 

S  +  N 

SB(T) 

S3(2T) 

S  B(T) 

S  +  N 

Note  that  since  the  matrix  (C+NI)  is  real  and  symmetric,  it  may 
be  factored  as  follows 


a11 

0 

0 

“  a11 

a21 

a31 

(C+NI)  =  a2. 

a22 

0 

0 

a22 

a32 

a31 

a32 

a33 

0 

0 

a33 
( A-7 ) 

where  the  unknown  coefficients,  a^ j ,are  all  real.  Defining 
the  matrix  factorization  as  (C+NI)  =  MM',  then  the  reader  may 
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easily  verify  the  following  solution. 
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11  “ 

^s  +  n 

S  T) 

11 

Vs  +  N 

S8(2T) 

31  “ 

<^S  +  N 

V(S  +  N)2  -  (  S  &(  T)  )  2 

22  * 

^S  +  N 

(SB(T)  [(S+N)  -  (  S  0(  2T 

32  * 

y/TTW 

*\f  S+N) 

3-[(S8(2T))2+2(S0(T))2; 

(A-8) 


V<S  V  n)2  -  ( S  6( T) ) 


The  sequence  of  complex  data  samples  {zq,  z^ ,  z2 }  is  then 
generated  using  the  following  matrix  linear  transformation 
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The  input  samples  gQ,  g , ,  g2  are  uncorrelated  complex 
Gaussian  random  variables 

90  “  *0  +  jy0 

9]  »  x,  +  jy1  (A-10) 

g2  ■  x2  +  jy2 

Thus,  to  obtain  the  properly  correlated  sequence  of  complex 
data  samples  we  generate  6  independent  (mean*0,  variance* 1/2) 
Gaussian  random  variables  xq,  yo,  xj,  yi ,  X2  and  y2 
and  obtain  the  corresponding  simulated  complex  radar  signal 
samples  as  indicated  above  in  equation  A-9.  It  is  easy  to  show 
that  the  complex  data  samples  {zQ,  z^ ,  *2 }  have  the  desired 

Hermitian  symmetric  covariance  matrix  (see  Equation  A-3). 

It  is  easy  to  verify  the  fact  that  the  variables  {zQ,  Zy  z2 } 

also  satisfy  the  required  circular  property  (see  reference  [6]) 

EU^Xj}  *  0  for  all  i,  j  (A-11) 

We  remark  that  this  property  will  not  be  satisfied  if  one  were 
to  take  the  input  data  vector  to  be  real  Gaussian  variates 
(e.g.,  gQ  *  xQ,  g1  =  x^  g?-  x2).  Finally,  it  is  clear  that 

the  method  described  in  this  appendix  generalizes  to  the 
general  L-pulse  case. 


APPENDIX  B  ERROR  ANALYSIS  OF  MEAN  VELOCITY  ESTIMATOR 
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This  appendix  provides  a  summary  of  the  error  analysis 
of  the  mean  velocity  estimator  for  burst  waveforms  which  are 
processed  using  the  Rummler  time  domain  approach.  Specifically 
we  have  considered  a  burst  waveform  comprised  of  L  contiguous 
pulse  returns,  denoted  {*0,  *2» * • * »*l-1 this 

contiguous  pulse  train  we  may  compute  estimates  of  the  mean 
doppler  frequency  using  various  lag  products.  For  example,  we 
might  construct  the  following  estimates. 


iy  «  -yir  ARG 

f2  ’  TlTTTf)  ARG 

A 

f3  "  i? nf)  ARG  1 


where 

r(t)  =  -x^r 
i(2T)  = 

i(3T)  » 


etc 


Lf 

i»0 

L-3 

I 

i=0 


z*i  zi+1 


zt  z. , _ 
i  i+2 


L-4 

I  z*  z. 

i*0 


i+3 


( B—  1  ) 


( B-2 ) 


etc. 
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Each  of  the  above  estimators  f^,  ^3/  etc  yields  valid 

estimates  of  the  mean  doppler  frequency.  One  is  interested  in 
determining  the  accuracy  of  these  estimates.  A  brief 
discussion  of  the  accuracy  analysis  is  given  in  the  following 
paragraphs.  This  analysis  is  based  upon  the  work  of  Reference 
[6]  . 

An  analysis  of  the  accuracy  of  the  Rummler  time  domain 
estimator  (see  appendix  A,  Reference  (6])  yields  the  following 

A 

error  expression  for  the  accuracy  of  the  estimator  f^ 

Var { ffc  }  = 


1 

—  - 

8^(KT)Z 


Re 


A 

R(KT)  |21  _ 

[/r(KT)  \2 

!ll 

R(KT)|  J  E 

Hr(KT)  J 

IJ 

(B-3) 


Omitting  the  lengthy  derivations  of  the  terms  of  the  above 
expression  (see  Reference  [6]  for  basic  approach)  one  may 
obtain  the  final  expression  shown  below. 

Var  {fR}  =  [l/8*2(KT)20(KT)  ] 

1  m=( L-K-1 )  ? 

— - 7  l  ( 0  (mT)- 0( (m+K)T) 0( (m-K)T) )  [L-K-lml] 

(L-K)  m=-( L-K-1 )  1  1 

+  TTTkT  £2<n/S)  +  (N/S)2  ] 


(B-4) 


(L-K) 

2 

(L-K) 


2  [(N/S)  0(  2KT)  (L-2K)  ] 

_ _ _ _ i 


L>2K,  otherwise=0 
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Note  that  the  third  term  in  the  above  expression  is  valid  only 
when  L>2K,  otherwise  this  term  is  zero. 

A  A 

Since  and  f 2  yield  two  viable  estimates  of  the 
mean  doppler  frequency,  one  may  attempt  to  improve  the  accuracy 
of  the  mean  doppler  frequency  estimate  by  taking  a  linear 

A  A 

combination  of  the  estimates  f1  and  f  2 . 

f 1 2  =  C1  f1  +  c2f2  (B"5) 
The  weights  (c-|,  C2)  could  be  selected  to  provide  a  better 
estimate  of  the  mean  doppler  frequency  (better  in  the  minimum 
variance  sense).  Selecting  uniform  weights  (c-|=C2BS1/2) 
might  also  be  considered  as  a  simpler  improved  estimator.  In 

A 

either  case,  the  accuracy  of  estimate  f y 2  is  given  by  the 
expression 


Var  U12}  =  c^  Var  {f  1  }  +  c\  Var  {f  2 } 

A  A 

+  2  c1  c2  Cov  {f 1 ,  f 2  } 

A 

Selecting  (c^ ,  c2 )  to  make  f12minimum  variance  and 
unbiased  we  obtain  the  optimal  solution 


( B-6 ) 


c*  =  var  tf-,1  -  Cov  tf2  ,  f,}  (B_7) 

Var  {f  1  }+Var  {f2  }-2  Cov{f2,  ^  } 

c*  =  1-c*  ( B-8 ) 

Similar  expressions  can  be  obtained  for  an  improved  estimator 
of  the  form 


r«jv  aJI \  W.r. 


^  A*.  Mj*  • 
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f  1 23  =  C1  f1  +  C2  f 2  +  C3  f3 


(B-9) 


To  evaluate  the  theoretical  accuracy  of  these  various 
estimators  it  is  clear  that  one  must  evaluate  the  covariance 


*  A  A 


between  the  errors  of  each  Rummler  estimate  f^,  fj,  f^, 
etc.  If  the  covariance  between  the  errors  of  these  estimates 
is  small  then  significant  improvement  in  mean  doppler  frequency 
estimates  can  be  achieved. 


The  covariance  between  estimate  fK  and  fj  (where  K>J ) 
is  evaluated  as  follows.  Extending  the  results  of  Appendix  A 
of  Reference  [6  ]  we  have  the  expression 


Cov  {fR,  fj}  =  [i/8*2{KT)(JT)  ] 


fJ  /r(JT)|  {R(KT)f  -  /r(JT)|  (R(KT) 
\R(W)/  \R[JT)J  \R(KT) 


where. 


R<JT>  - 


R( KT }  =  1 

L-K 


L-J-1 


L-K- 1 


2i  Zi+J 


iio 


B- 1  0 ) 


(B-11) 


Omitting  the  lengthy  derivations,  the  result  may  be  shown  to  be 
given  by  the  following  expression. 
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A  A  a 

COV  *fK'  fjJ  =  [V8it2(KTMJT)(L-J)(L-K)  0(KT)  3(JT)  ] 

L-J-1  L-K-1 

i~0  i-0  e((  j~i)T)  0((i-j+J-K)T)-0((j-i+K)T)  0((i-j+J)T). 

( B- 1 2 ) 

+  2 ( L-K )  0( ( J-K)T) (N/S) 

-  2( L-K-J ) 0( ( J+K ) T) ( N/S ) 

L-K-J  >  0,  otherwise  =  0 


Note  that  the  last  term  in  the  above  expression  is  valid  only 
when  L-K-J  >  0f  otherwise  this  term  is  zero. 


APPENDIX  C  ERROR  ANALYSIS  OF  SPREAD  ESTIMATOR 

This  appendix  provides  a  summary  of  the  error  analysis 
equations  used  to  evaluate  the  performance  accuracy  of  the 
Rummler  spread  estimator.  The  basic  algorithm  is  defined  by 
the  following  expression 


W  =  — - —  1  -  I  p ( T )  I  1/2  (C-1) 

T  '  1 

a 

where  p(T)  is  the  estimated  correlation  coefficient  of  the 
underlying  signal  process.  The  correlation  coefficient  is 

estimated  from  the  sample  correlation  function  as  follows 


P(T)  =  — —  ( C-2 ) 

A 

R( 0 )  -  N 

Thus,  it  is  assumed  that  the  noise  power,  N,  is  known  or  can  be 
accurately  estimated  from  range  gate  data  comprised  of 
noise-only  samples.  Also,  we  have  the  terms 


( C-3 ) 


*<T>  “  ^  zi+1  <c-4» 

This  estimator  provides  reasonably  good  estimates  of  the  spread 
parameter  "W"  for  small  values  of  T.  To  see  this,  one  simply 
considers  for  a  large  number  of  pulses  (L+«) 
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R(T)  ^ 
Thus,  for  large  L, 


S  3(T)e 


p(T)  e 


i2  *fT 


-2if2W2T2 


(C-6) 


(C-7 ) 


Approximating  the  right  hand  side  of  the  above  equation  by  the 
term  1-2*2W2T2  we  obtain. 


y[2*T 


2  2  2 
2«*VTT* 


(C-8 ) 


Thus,  for  small  T  we  have 


This  derivation  was  not  intended  to  be  rigorous  but  was 
presented  simply  to  demonstrate  that  the  algorithm  W  may 
provide  reasonably  good  estimates  of  the  spread  parameter,  W, 
for  small  T.  This  algorithm  is  however,  a  biased  estimator 
of  the  parameter  W  since 

e{w}  *  W  (C-10) 

and  thus  the  analysis  of  the  performance  accuracy  of  this 

algorithm  must  include  a  determination  of  both  its  bias  and 

variance  terms.  It  may  be  shown  (see  Reference  [10]  for 

details  of  the  derivation)  that  the  variance  and  bias  errors 

for  the  Rummler  spread  estimator  are 


^  -  A  tv*-.  v  Au.-  v.-.i  v  v  o  ^  o(V 
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A.  Random  Error  Equation 
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Bias  Error  Equation 
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+  0.25  ({TERM  3}  -  {TERM  4}) 
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where 


y[2*T 


1-B(T) 


(C- 1 7 ) 


Another  algorithm  for  estimating  the  spread  parameter, 
W,  is  the  ’  jarithmic  version  of  Rummler's  algorithm,  defined 
by  the  following  expression  (see  Reference  [10]). 


W.  =  — • - —  I-  In  |  P(T) 

109 


(C-18) 


This  algorithm  is  reported  to  have  less  bias  than  the  Rummler 


version,  especia 


for  the  case  where  the  signal  process  has  a 


power  spectral  uensity  which  is  Gaussian  in  shape  (and  hence,  a 
Gaussian  autocorrelation  function).  Taking  the  same  approach 
as  for  the  previou  iread  estimator  it  is  easy  to  determine 
the  error  equatons  for  the  logarithmic  spread  estimator. 


:  •  - . -ji  ■ 


The  variance  and  bias  errors  for  this  estimator  are  summarized 
below 


where 


(-lntf(T))/2 
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Fig.  1.  Normalized  standard  deviation  of  mean  velocity  estima¬ 
tors  versus  normalized  pulse  spacing  (16  pulses  total,  SNR  =  -6 
dB) . 
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BURST 


CR  BOUND 


Standard  deviation  (m/s)  of  mean  velocity  estimators 
V3  versus  pulse  spacing  (u  sec)  (16  pulse  purst,  SNR 
<*v  *  150  m/s,  f  -  9.1  GHz) . 
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Fig. ^6.  Standard  Deviation  (m/s)  of  mean  velocity  estimators 

71'  V12'  V1 23  versus  Pulse  spacing  (y  sec)  with  optimal  weights 
(16  pulse  burst,  SNR  =  -6  dB,  =  150  m/s,  f  =  9.1  GHz). 


STANDARD  DEVIATION  (m/s) 


STANC 


Fig.  8.  Standard  Deviation  (m/s)  of  velocity  spread  estimators 
versus  signal-to-noise  ratio  (dB)  (16  pulse  burst,  150  ir./s, 

f  =  9.1  GHz,  T  =  5  usee) . 
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algorithm  for  estimating  the  spectral  mean  is  derived.  This  M.L.E  algorithm  is  shown  to  include,  as  a  special  case, 
the  spectral  mean  estimator  originally  proposed  by  R.W.  Miller  [1]  for  processing  of  burst  waveforms.  Also,  when  the 
burst  waveform  is  a  simple  pulse  pair,  the  M.L.E.  algorithm  reduces  o  'he  spectral  mean  estimator  originally  pro¬ 
posed  by  W.D.  Rummler  (2).  The  Cramer-Rao  bound  for  estimating  the  ..pectral  mean  using  burst  waveforms  is  also 
derived.  Simplifications  to  the  exact  maximum-likelihood  algorithm  are  proposed  and  the  performance  of  various 
estimators  is  compared  to  the  Cramer-Rao  lower  bound.  Some  preliminary  results  of  studies  of  spectral  width  estima¬ 
tors  are  also  presented. 
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